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Abstract 

The existence of a ne'w force beyond the Standard Model is compelling because it could explain several striking astrophysical 
observations 'which fail standard interpretations. We searched for the light vector mediator of this dark force, the U boson, 'with 
the KLOE detector at the DAONE e^e" collider. Using an integrated luminosity of 1.54 fb *, 'we studied the process e*^e“ —» Uy, 
■with U —> e^e , using radiative return to search for a resonant peak in the dielectron invariant-mass distribution. We did not find 
evidence for a signal, and set a 90% CL upper limit on the mixing strength bet'ween the Standard Model photon and the dark photon, 
at 10“®-10“"* in the 5-520 MeV/c^ mass range. 
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1. Introduction 

The Standard Model (SM) of particle physics has received 
further confirmation 'with the discovery of the Higgs boson [lH- 
13, ho'wever, there are strong hints of physics it cannot explain, 
such as neutrino oscillations ei and the measured anomalous 
magnetic moment of the muon 0. Eurthermore, the SM does 
not provide a dark matter (DM) candidate usually advocated 
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as an explanation of the numerous gravitational anomalies ob¬ 
served in the universe. Many extensions of the SM ||7j? - [TOl 
consider a Weakly Interacting Massive Particle (WIMP) as a 
viable DM candidate and assume that WIMPs are charged un¬ 
der a ne'w kind of interaction. The mediator of the ne'w force 
■would be a gauge vector boson, the U boson, also referred to 
as a dark photon or A'. It 'would be produced during WIMP 
annihilations, have a mass less than t'wo proton masses, and a 
leptonic decay channel in order to explain the astrophysical ob¬ 
servations recently reported by many experiments lfTTH2T1l . 

In the minimal theoretical model, the U boson is the light¬ 
est particle of the dark sector and can couple to the ordinary 
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SM photon only through loops of heavy dark particles charged 
under both SM U(1 )y and dark U(1 )d symmetries ||22j? - l26l . 
These higher-order interactions would open a so-called kinetic 
mixing portal described in the theory by the Lagrangian term 
Xmix = - f ^,7''^^Dark’ ^here Ff^ is the SM hypercharge gauge 
field tensor and is the dark field tensor. The e parameter 
represents the mixing strength and is the ratio of the dark and 
electromagnetic coupling constants. In principle, the dark pho¬ 
ton could be produced in any process in which a virtual or real 
photon is involved but the rate is suppressed due to the very 
small coupling (e < 10“^). In this respect, high-luminosity 
(9(GeV)-energy e^e“ colliders play a crucial role in dark pho¬ 
ton searches Il27l429l . 

We investigated the e^e“ —> Uy process by considering the 
U boson decaying into e^e". At the level of coupling accessi¬ 
ble by KLOE in this channel the U boson is expected to decay 
promplty leaving its signal as a resonant peak in the invariant- 
mass distribution of the lepton pair. The energy scan was per¬ 
formed by applying the radiative return method which consists 
of selecting the events in which either electron or positron emits 
an initial-state radiation (ISR) photon which carries away a part 
of the energy and effectively changes the amount of the en¬ 
ergy available for U boson production. The selected initial- 
and final-state particles are the same as in the radiative Bhabha 
scattering process so we receive contributions from resonant s- 
channel, non-resonant f-channel U boson exchanges, and from 
s-t interference. The finite-width effects related to i-channel 
annihilation sub-processes, scattering f-channel and s-t inter¬ 
ference are of order of Tu/my for the integrated cross section 
and can be neglected with respect to any potential resonance 
we would observe; Tu ~ 10“^-10“^ MeV for the coupling 
strengths to which we are sensitive lIMI . The non-resonant f- 
channel effects would not produce a peak in the invariant-mass 
distribution but could, in principle, appear in analyses of angu¬ 
lar distributions or asymmetries. We are going to report exclu¬ 
sively on resonant s-channel U boson production. 

Using a sample of KLOE data collected during 2004-2005, 
corresponding to an integrated luminosity of 1.54 fb *, we de¬ 
rived a new limit on the kinetic mixing parameter, e^, approach¬ 
ing the dielectron mass threshold. 


read out at both ends by 4880 photomultiplier tubes, ultimately 
resulting in an energy resolution of cteIE - 5.1%I -\jE(GeV) 
and a time resolution of cr, = 57 ps/ -\jE(GeV) © 100 ps. A su¬ 
perconducting coil around the EMC provides a 0.52 T field to 
measure the momentum of charged particles. A cross sectional 
diagram of the KLOE detector is shown in Eigure[T] 

The trigger lOTI uses energy deposition in the calorimeter 
and drift chamber hit multiplicity. To minimize backgrounds 
the trigger system includes a second-level cosmic-ray muon 
veto based on energy deposition in the outermost layers of the 
calorimeter, followed by a software background filter based on 
the topology and multiplicity of energy clusters and drift cham¬ 
ber hits to reduce beam background. A downscaled sample is 
retained to evaluate the filter efficiency. 



Figure 1: Cross section of the KLOE detector. 


2. KLOE detector 

The Erascati cf) factory, DA<I)NE, is an e^e“ collider running 
mainly at a center-of-mass energy of 1.0195 GeV, the mass of 
the 4) meson. Equal energy electron and positron beams collide 
at an angle of ~25 mrad, producing cf) mesons nearly at rest. 

The KLOE detector consists of a large cylindrical Drift 
Chamber (DC) ^ with a 25 cm internal radius, 2 m outer 
radius, and 3.3 m length, comprising ~56,000 wires for a total 
of about 12,000 drift cells. It is filled with a low-Z (90% he¬ 
lium, 10% isobutane) gas mixture and provides a momentum 
resolution of crp^lp^_ ss 0.4%. The DC is surrounded by a 
lead-scintillating fiber electromagnetic calorimeter (EMC) If32l 
composed of a cylindrical barrel and two end-caps providing 
98% coverage of the total solid angle. Calorimeter modules are 


3. Event selection 

Using 1.54 fb ' of KLOE data we have searched for U bo¬ 
son production in the process e^e“ ^ Uy followed by U ^ 
e^e . The center-of-mass energy of the collision depends on 
the amount of energy carried away by the initial-state radia¬ 
tion (ISR) photon. The irreducible background originates from 
the e^e“ —> e^e'y radiative Bhabha scattering process, having 
the same three final-state particles. The reducible backgrounds 
consist of e^e“ ^ p-^p^y, e^e“ —» Tt^Tt^y, e^e" —» yy (where 
one photon converts into an e^e“ pair), and e^e“ ^ cf) —» 
pTT*^ —> Tt^TT^TT^, as Well as other cf) decays. The expected U bo¬ 
son signal would appear as a resonant peak in the invariant- 
mass distribution of the e^e“ pair, mee. This search differs from 
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the previous KLOE searches Il?4l43^ in its capability to probe 
the low mass region close to the dielectron mass threshold. 

We selected events with three separate calorimeter energy de¬ 
posits corresponding to two oppositely-charged lepton tracks 
and a photon. The final-state electron, positron, and photon 
were required to be emitted at large angle (55° < 6 < 125°) 
with respect to the beam axis, such that they are explicitly de¬ 
tected in the barrel of the calorimeter, see Figure [T] The large- 
angle selection greatly suppresses the f-channel contribution 
from the irreducible Bhabha-scattering background which is 
strongly peaked at small angle. Since we are interested mostly 
in the low invariant-mass region, we select only events with a 
hard photon, Ey > 305 MeV, chosen to select a subsample of 
the events generated by our MC simulation. We required both 
lepton tracks to have a first DC hit within a radius of 50 cm 
from the beam axis and a point-of-closest-approach (PCA) to 
the beam axis within the fiducial cylinder, < 1 cm and 
-6 < < 6 cm, entirely contained within the vacuum pipe 

eliminating background events from photons converting on the 
vacuum wall. We eliminated tightly spiralling tracks by re¬ 
quiring either a large transverse or a large longitudinal mo¬ 
mentum for each of the lepton tracks, > 160 MeV/c or 
Pj. > 90 MeV/c. We require that the total momentum of the 
charged tracks is (|pe+| + IPe l) > 150 MeV/c to avoid the pres¬ 
ence of poorly reconstructed tracks. A pseudo-likelihood dis¬ 
criminant was used to separate electrons from muons and pi- 
ons El. A further discrimination from muons and pions was 
achieved using the Mtrack variable. Mtrack is the X mass for an 
X^X~y final state, computed using energy and momentum con¬ 
servation, assuming mx+ = mx- llJTl . In Figure]^ the Mtrack 
distribution is reported for measured data and for all the rele¬ 
vant MC simulated background components. Including the cut 
Mtrack < 70 MeV/c^ we were left with 681,196 events at the end 
of the full analysis chain. 



^track (MeV/c^) 


Figure 2: (color online) The track mass distribution before event selection for 
measurement and expected background simulations. Some of the simulations 
had a prescaling for Mtrack < 80 MeV/c^, which has been accounted for in the 
background evaluation. The measurement data set was prescaled for Mtrack > 
85 MeV/c^. The track mass variable peaks at the mass of the charged track in 
the final state for events with two charged tracks and a photon. The selection 
region is Mtrack < 70 MeV/c^. 


4. Simulation and efficiencies 

We used MC event generators interfaced with the full KFOE 
detector simulation, GEANFI including detector resolu¬ 
tions and beam conditions on a run-by-run basis, to estimate 
the level of background contamination due to all of the pro¬ 
cesses listed in the previous section. Excluding the irreducible 
background from radiative Bhabha scattering events, the con¬ 
tamination from the sum of residual backgrounds after all anal¬ 
ysis cuts is less than 1.5% in the whole range, and none 
of the background shapes are peaked, eliminating the possibil¬ 
ity of a background mimicking the resonant U boson signal. 
The irreducible Bhabha scattering background was simulated 
using the Babayaga-NLO Il39ti42l event generator implemented 
within GEANFI (including the s-, t-, and s-t interference chan¬ 
nels) and is shown in Figure along with the measured data 
after subtracting the non-irreducible background components. 
No signal peak is observed. 



Figure 3: (color online) Dielectron invariant-mass distribution from mea¬ 
surement data with non-irreducible backgrounds subtracted compared to the 
Babayaga-NLO MC simulation. 

In order to evaluate the U boson selection efficiency we used 
a modified version of the Babayaga-NLO event generator imple¬ 
mented within GEANFI, such that the radiative Bhabha scat¬ 
tering process was only allowed to proceed via the annihila¬ 
tion channel, in which the U boson resonance would occur. In 
order to create a large-statistics sample in our region of inter¬ 
est we restricted the Babayaga-NLO generated events to within 
50° < < 130° and > 300 MeV. The generator- 

level efficiency due to this restriction was evaluated using a 
Phokhara MC simulation ll43ll . The total efficiency is eval¬ 
uated as the product of the generator-level efficiency and the 
event-selection efficiency, containing the cuts in Section|^con- 
ditioned to the generator-level restriction as well as the trigger 
efficiency, and is shown in Figure The decrease in efficiency 
as mee ^ 2me comes from the requirement on the total momen¬ 
tum of the charged tracks. 


3 













0.04 


cross section, 



Figure 4: Smoothed distribution of the total efficiency defined as the product 
of the selection efficiency for the e‘''e“ —> Uy —> e‘''e“y final state evaluated 
using the Babayaga-NLO event generator modified to allow only the ^-channel 
process, and the generator-level efficiency evaluated from a Phokhara MC sim¬ 
ulation. 

5. Upper limit evaluation 

We used the CLs technique ll4^ to determine the limit on the 
number of signal U boson events, Nv, at 90% confidence level 
using the mee distribution. The invariant-mass resolution, 
is in the range 1.4 < < 1.7 MeV/c^. Chebyshev poly¬ 

nomials were fit to the measured data (+15crm^J, excluding the 
signal region of interest (+3cr„_,J. The polynomial with;^f^/Adof 
closest to 1.0 was used as the background. A Breit-Wigner peak 
with a width of 1 keV smeared with the invariant-mass resolu¬ 
tion was used as the signal. An example of one specific CLg 
result is shown in Figure]^ yielding an upper limit of N\j = 215 
U boson events at m\j - 155.25 MeV/c^ at the 90% confidence 
level. The x^l^dof was 1.09 for this Chebyshev-polynomial 
sideband fit. 



UL[cr(e^e ^ Uy,U ^ e^e )] = ^ (i) 

L fett 

where L is the luminosity and egff is the total selection effi¬ 
ciency. The limit is shown in Figure]^ 



Figure 6: Upper limit on the cross section crfe^e —» Uy, U -+ e^e ). 

We then translated the limit on Nd to a 90% confidence level 
limit on the kinetic mixing parameter as a function of mge as 
in I36l, 


2 , ^ Au(OTeg) _1_ 

eeff(?Mee) I(mee) L 


( 2 ) 


where the radiator function Hiniee) was extracted from 

dcrgg.y/dmee = H{mee, s,cos(0-),)) ■ cr2;™(mge) 

using the Phokhara MC simulation Il43l to determine the ra¬ 
diative differential cross section, /(mgg) is the integral of the 
cross section (T(e^e“ —> U ^ e^e^), L = 1.54 fb“* is the inte¬ 
grated luminosity, and egfj (mgg) is the total efficiency described 
in Section]^ Our limit is shown in Figure [7] along with the in¬ 
direct limits from the measurements of (g - 2)g and (g - 2)^ at 
5cr shown with dashed curves. Limits from the following di¬ 
rect searches are shown with shaded regions and solid curves: 
E141 HD, E774 Ha, KLOE(4) ^ pU, U ^ e+e-) iMlEa, 
Apex ESI, WASA E3, HADES ESI, A1 ESI, KLOE(e+e- ^ 
Uy, U ^ y p.-) ESI, BaBar ESI, and NA48/2 BTl . 


Figure 5: (color online) The CLs result at 90% CL for mu = 155.25 MeV/c^ 
showing the measured data, the Chebyshev-polynomial sideband fit, and the 
signal shape scaled to the CLs result. 


The upper limit at 90% confidence level on the number of 
U boson events, UL(Au), can be translated into a limit on the 
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Figure 7: (color online) Exclusion limits on the kinetic mixing parameter 
squared, s^, as a function of the U boson mass. The red curve labeled KLOE( 3 ) 
is the result of this article while the curves labeled KLOE(i) and KLOE( 2 ) indi¬ 
cate the previous KLOE results. Also shown ai‘e the exclusion limits provided 
by E141, E774, Apex, WASA, HADES, Al, BaBar, and NA48/2. The gray 
band delimited by the dashed white lines indicates the mixing level and mu 
parameter space that could explain the discrepancy observed between the mea¬ 
surement and SM calculation of the muon (g-2),^. 

6. Systematic uncertainties 

The background was determined by Chebyshev-polynomial 
sideband fits. The parameters of the polynomials were then var¬ 
ied within Icr to determine the maximum variation of the poly¬ 
nomial shape. The uncertainty of each bin was set to the extent 
of that variation evaluated at the bin center. An example of 
the error bars on the Chebyshev-polynomial sideband fits can 
be seen in Figure]^ These bin uncertainties were taken into ac¬ 
count in the CLs procedure when determining VcLs(^ee)- Since 
the irreducible background is smooth for each fit range, we 
assume the Chebyshev polynomials sufficiently represent the 
background with negligible systematic uncertainty. Any uncer¬ 
tainty in the shape of the smeared resonant peak was also taken 
to be negligible. 

The efficiency of the e^e“ ^ e^e“Y event selection was 
determined by taking the ratio of the set of simulated events 
that passed the selection criteria to the total simulated sample. 
We apply a 0.1% systematic uncertainty due to the Babayaga- 
NLO event generator ll39H42l . a 0.1% systematic uncertainty 
for the trigger, and a 0.1% systematic uncertainty for the soft¬ 
ware background filter. All together the uncertainty on the se¬ 
lection efficiency is dominated by the statistical uncertainty on 
the selected sample. A Phokhara MC simulation P3ll was per¬ 
formed to evaluate the generator-level efficiency due to the re¬ 
striction > 300 MeV and 50° < < 130°. The selec- 

y e+,e 

tion efficiency and the generator-level efficiency are combined 
to give the total efficiency, eeff(?Mee)- The uncertainty is given 
as the error band in Figure]^ again dominated by the statistical 
uncertainties in the simulated data set. 

There are two effects that contribute to the uncertainty in the 
radiator function, //(mee). First, since the value of is 

taken from simulated data, we must take into account the statis¬ 


tical uncertainty on those values. Second, we assume a uniform 
0.5% systematic uncertainty in the calculation of //(mee), as 
quoted in Il43l I52l - l54l . The uncertainty in the integrated lumi¬ 
nosity is 0.3% El. The uncertainties on eefr('Wee). and 

L, propagate to the systematic uncertainty on eP'imee) via 0. A 
summary of systematic uncertainties is presented in Tabled] 

Table 1: Summary of systematic uncertainties. The uncertainties on the effi¬ 
ciency, radiator function, and cross-section integral vary as a function of 
The numbers quoted here correspond to the largest estimate within our mge 
range. 

Systematic source Relative uncertainty 

Background (sideband fit) negl. 

fefF (mee) 2% 

MC generator, 0.1% 

Trigger, 0.1% 

Software background filter, 0.1% 

Event selection, 2% 

//(mee) 0.5% 

/(mee) negl. 

L 0.3% 


7. Conclusions 

We performed a search for a dark gauge U boson in the pro¬ 
cess e^e“ ^ Uy with U —> e^e“ using the radiative return 
method and 1.54 fb ' of KLOE data collected in 2004-2005. 
We found no evidence for a U boson resonant peak and set a 
90% CL upper limit on the kinetic mixing parameter, e^, at 
10“®-10“"^ in the U-boson mass range 5-520 MeV/c^. This 
limit partly excludes some of the remaining parameter space 
in the low dielectron mass region allowed by the discrepancy 
between the observed and predicted (g-2)^^. 
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